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We show how effectively effective quantum field theories work in nuclear physics. Us- 



ing the physically transparent cut-off regularization, we study the simplest nuclear systems 
of two nucleons for both bound and scattering states at a momentum scale much less than 
the pion mass. We consider all the static properties of the deuteron, the two-nucleon scat- 
tering phase-shifts, the n + p — > d + ^f process at thermal energy and the solar proton fusion 
process p + p —* d + e + + v e , and we demonstrate that these are all described with great 
accuracy in the expansion to the next-to- leading order. We explore how a "new" degree 
of freedom enters in an effective theory by turning on and off the role of the pion in the 
Lagrangian. 
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1 Introduction 



The idea of effective quantum field theory]!], |2], |3| is extremely simple. At low energy 
where non-perturbative effects of QCD dominate, the relevant degrees of freedom are not 
quarks and gluons but hadrons. Let the relevant hadronic degrees of freedom be represented 
by the generic field <£. Separate the field into two parts; the low-energy part $l in which 
we are interested and the high-energy part &h which does not interest us directly. We 
delineate the two parts at the momentum scale characterized by a cutoff, say, A]_. In the 
generating functional Z that we want to compute, integrate out the <&h field and define an 
effective action S £ ^[^l\ = i In ( f{d&H]e lS ^ L '^ H ^) • Then what we need to compute is 




Given the right degrees of freedom effective below the cutoff Ax, we expand the effective 
action as S^f = J2T where Qjf^iJ are local field operators allowed by 

the symmetries of the problem, and Cj's are constants, which are required to satisfy the 
"naturalness condition". The Qj[^i]' s are ordered in such a manner that, for low-energy 
processes, the importance of the term diminishes as i increases. The contribution 

of the <&h degrees of freedom integrated out from the action is not simply discarded but 
instead gets lodged in the coefficients Cj as well as in higher-dimensional operators. The 
strategy of effective field theories is to truncate the series at a manageable finite order, and 
obtain desired accuracy to capture the essence of the physics involved. The separation of 
<J> into and &h involves two steps: The first is to "decimate" the degrees of freedom, 
that is, to choose the fields relevant for the physics to be incorporated explicitly into 
the Lagrangian, and the second is a regularization procedure which regulates the high- 
momentum contribution. Let m be the lightest mass of the degrees of freedom that are 
integrated out. Clearly the most physically transparent regularization is to introduce a 
momentum cutoff with Ai ~ m. A different cutoff introduces different coefficients C{ for 
the given set of fields. The idea is then to pick a value of the cutoff that is low enough to 
avoid unnecessary complications but high enough to avoid throwing away relevant physics. 

In the full (un-truncated) theory, the location of the cutoff is entirely arbitrary, so 
we could have chosen any values of the cutoff other than Ai without changing physics. The 
cut-off dependence of Qj^l] is compensated by the corresponding changes in d through 
the renormalization group equation. In such a theory, which is rarely available in four 
dimensions, physical observables should be strictly independent of where the cutoff is set. 
This of course will not be the case when the series is truncated. However, for a suitably 
truncated effective theory to be predictive, the observables should be more or less insensitive 
to where the cutoff is put. Should there be strong cutoff dependence, it would signal that 
there is something amiss in the scheme. For instance, it could be a signal for the presence 
of new degrees of freedom or of "new physics" that must be incorporated and/or for the 
necessity of including higher-order terms. How "new physics" enters in our case will be 
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studied in terms of the role that the pion - a "new degree of freedom" - plays at very 
low energy. Here we of course know exactly what the "new physics" is but there will be 
a lesson to learn from this exercise in exploring beyond "standard models" . The following 
"toy model" illustrates our point. Let us imagine that we knew that QCD is the correct 
microscopic theory for the strong interactions with chiral symmetry realized in the Nambu- 
Goldstone mode but the pions were not yet discovered and that the only "effective" degrees 
of freedom we had at hand were the proton and the neutron. Now the question to ask is 
how to set up a systematic theory that could describe very low-energy properties of the 
deuteron, the NN scattering and the electroweak matrix elements and in doing so how to 
infer that a "new degree of freedom" (pion) lurked at ~ 140 MeV. We will see in Section 6 
how this question could be studied in the real low-energy world of nuclei. 

An important aspect of effective field theory is that the result should not depend 
upon what sort of regularization one uses. In principle, therefore, one could equally well 
use the dimensional regularization or the cut-off regularization or something else as long as 
one makes sure that all relevant symmetries are properly implemented W^. Calculated to 
high enough order in a consistent expansion, the results should in principle be the same up 
to that order although in practice there is a certain amount of art involved in the choice 
of the scheme. This issue in the context of the dimensional regularization and its variants 
was addressed recently in Q. Here we will use the cut-off regularization for the reason 
that the procedure is simple and physically transparent |5|, |6| when limited to the order we 
shall adopt, namely, to the next-to-leading order (NLO). For some of the problems treated 
here such as neutron-proton (np) scattering, one could also use in a predictive manner a 
modified dimensional regularization in which certain power divergences are removed (i.e., 
PDS scheme) J?]] . For further comments on the PDS scheme, see Gegelia It is not yet 
fully clear to us how such a scheme despite its advantages could avoid some of the problems 
raised in f§ |, §. 

In this paper, we show that an effective field theory with finite cut-off regulariza- 
tion works surprisingly well for all low-energy observables in the two-nucleon bound and 
scattering states; e.g., the static properties of the deuteron, the low-energy np scatter- 
ing, the radiative np capture, n + p — > d + 7, at thermal energy and the proton fusion, 
p + p — ► d + e + + v e , in the Sun. Some of the results were reported in previous pub- 
lications m & [K]] but they will be included here both for the sake of completeness and 
for presenting our more recent sharpened thoughts on the matter. One of the significant 
results of this paper is a quantitative justification from "first principles" of the previously 
formulated "hybrid method" |8| for the radiative np capture process n + p — > d + 7. We 
justify quantitatively the procedure of using the high-quality phenomenological two-nucleon 

* lr rhere are of course certain issues which raise technical, though not fundamental, problems, e.g., how to 
preserve chiral symmetry in the cut-off procedure and how to correctly handle short-distance physics in the 
standard dimensional regularization etc. For a recent discussion on this matter, see M. 
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wavefunctions ^ supplied by the Argonne vis potential O] in calculating the matrix ele- 



ment of the electroweak operators derived by chiral perturbation theory with "irreducible" 



diagrams, a procedure first suggested by Weinberg [12| (for an early implementation, see 
[HH). Similar support is given for the hybrid- method calculation based on the high-quality 
phenomenological wavefunctions of the proton fusion process p + p^d + e + + v e crucial in 
the solar neutrino problem. This constitutes an important application of effective theories 
for nuclei to astrophysics. 

2 Formulation of the Problem 

The problem we shall address is the following. Suppose we wish to calculate, starting 
from the "fundamental principles" of QCD, the structure of the deuteron and the scattering 
of two nucleons at very low energy, with the characteristic momentum probed p much less 
than the pion mass ~ 140 MeV. What are the relevant variables? Quarks and gluons are 
clearly not the right variables. We introduce instead the nucleon field ip in the Lagrangian 
as a matter field Next we have the option of having or not having the pion field figure 



explicitly in the Lagrangian. In the real world, the pion has a mass of about 140 MeV, 
much larger than the momentum scale we are interested in but in the chiral limit, it would 
be degenerate with the vacuum. So one might ask whether it is necessary to retain the 
pion in the framework even when the pion mass is taken into account. To address this 
question, we will work to the next-to-leading order (NLO) in the momentum expansion a 
la Weinberg |12||.^ It is then sufficient to take the two- nucleon potential written to NLO 



and insert it into an infinite series of reducible graphs to be summed - which is equivalent 
to solving the Schrodinger or Lippman-Schwinger equation with that potential. 

For the sake of completeness of presentation, we first describe the general features of 
the counting rules adopted in this article. The way these counting rules are implemented 
in our formulation will become clear as we describe how the terms that enter into the 
observables are computed. 

For processes whose typical momentum scale can be characterized by the pion mass, 
the pion field should figure explicitly in the Lagrangian. With nucleons treated as "heavy" 



* 2 The latest phase-shift analyses of the two-nucleon data and the construction of phenomenological po- 
tentials to reproduce the resulting phase shifts have reached such a level that all the existing data below 
~350 MeV can be reproduced with high precision JIt], C3. We refer to these state-of-the art nucleon-nucleon 
potentials as high-quality phenomenological potentials. The two-nucleon wave functions generated by these 
potentials are referred to as high-quality phenomenological wave functions. 

#3 In principle we could start with a Lagrangian consisting of the Goldstone meson fields only, obtain the 
baryons as solitons ( "skyrmions" ) and then proceed to compute the fluctuations. But this line of research 
has not yet met with a quantitative success. See jl4| for a discussion on this matter. 

#4 Recently, Weinberg's counting rule (^) has been questioned by the authors of j?J in connection with the 
large scattering length for the S-wave np scattering. Our use of the Weinberg counting is not affected by 
this criticism. This matter is discussed in Appendix. 
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and the anti-nucleon degrees of freedom integrated out, one obtains a consistent and sys- 
tematic expansion scheme - heavy-baryon chiral perturbation theory - where nucleons and 
pions are the dynamical fields. Let Q be the typical momentum scale of the process or the 
pion mass, which is regarded as "small" compared to the chiral scale Ay ~ 1 GeV. The 



counting rule given by Weinberg! 12 1 is that, a Feynman graph comprised of A nucleons, Ne 
external (electroweak) fields, L irreducible loops and C-separated pieces is order of 0(Q U ) 
with 

ChPT: v = 2L + 2(C -l) + 2- (A + N E ) + Y,"i, (2) 

i 

n = di + y + a - 2 (3) 

where we have characterized each vertex i by d, the number of derivatives and/or the power 
of the pion mass, n« that of nucleon fields and e% that of external (electroweak) fields. The 



quantities v\ and C are defined so that i>i > even in the presence of external fields [15] 
and (C — 1) = for connected diagrams. 

Now if we limit ourselves to a very low- momentum regime, say, p <C f%, then even 
the pion field can be integrated out and we are left only with the nucleon field. Since we 
have integrated out the anti-nucleon degrees of freedom as well as all other meson fields, 
there are no irreducible loops and the corresponding counting rule can be obtained simply 
by putting the number of loops L equal to zero, 

EFT: v = 2{C-l) + 2-{A + N E )+Y J Vi (4) 

i 

with the same definition for Vi given in eq.(^). But the meaning of d\ in eq.(^) is changed: 
Since the pion field is integrated out, the di stands here only for the number of derivatives, 
and not the power of the pion mass. Then, up to NLO, the resulting potential consists of 
contact interactions and their two spatial derivatives. 

Electromagnetic interactions have different counting rules and should be treated sep- 
arately. Since the Coulomb interaction (between protons) is given as -% (where a ~ 1/137 
is the fine-structure constant and q is the momentum transferred), it is of relevance in an 
extremely small momentum region while it becomes irrelevant in other regions due to the 
smallness of a. Since we are interested in the proton fusion at threshold, we will explicitly 
include the Coulomb potential in the pp sector. 

If we increase the typical probe momentum to a value close to but lower than the 
pion mass, the pion degree of freedom will become "marginal" and its explicit presence will 
improve the convergence of the series. In this article, we will perform an NLO calculation 
both with and without the pion field, thereby revealing its interplay in the effective field 
theory (EFT). To this order, pion loops do not enter and hence we can simply focus on the 
(irreducible) potential. Our unregulated potential V has the form 



g\ o\ ■ qa 2 • q 47r - 



4/2 q2+m2 M 



V(q) = -n ■ r 2 |^ \" V + T7 C o + {C 2 P> + D 2 aV)qW + Z X Z 2 ^ (5) 



o 



q 
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with 

v M±^-fU., 2 ), ( 6) 

where q is the momentum transferred, qa — 1-26 the axial-vector coupling constant, f n ~ 
93 MeV the pion decay constant, M ~ 940 MeV the nucleon mass, and Z\Z% = 1 for the 
pp channel and zero otherwise. For the proton fusion process, we will also consider 0(a 2 ) 
corrections, i.e., the vacuum polarization (VP) potential and the two-photon-exchange (C2) 
potential; as in |1C], these will be treated perturbatively. The parameters C"s and D2 are 
defined for each channel. The D2 term, however, is effective only for the spin-triplet channel. 
Thus, there are two parameters (Co and C2) for each of the pp and np 1 So channels, and 
three (Co, C2 and D2) for the 3 <Si channel. 

It is perhaps useful to explain how the potential (||) is to be understood. The first 
(nonlocal) term is the pion exchange involving the Goldstone boson and hence completely 
known. The (local) terms in the square brackets represent the effect inherited from the 
degrees of freedom that have been integrated out, for instance the heavy mesons whose 
masses are much higher than the cutoff. It should be noted that when the pion is also 
integrated out, its effect will appear in the local terms, modifying the constants. This 
means that even if we drop the first term of (||), we still have part, if not most, of pionic 
effects in the theory. 

Inserted into an infinite bubble sum required to describe bound and quasi-bound 
states, the potential (||) will generate an infinite series of divergent terms unless the integrals 
are suitably regularized. This is expected in an effective theory which is not renormalizable 
in the conventional (Dyson) sense. As announced we shall regularize the divergences with 
a momentum cutoff. In || where only nucleonic contact terms and their derivatives were 
taken into consideration, a regularization that preserves the separability of the potential 
was used. Of course, the precise form of the regulator should not matter ||. Since both 
the Coulomb interaction and the one-pion-exchange potential are non-separable, we shall 
adopt a slightly different form of the regularization: in coordinate space, 

Viv) |^ e ^ r S A (q 2 )V(q), (7) 

where 5a (q 2 ) is the regulator with a cutoff A. For our purpose it is convenient to take the 
Gaussian regulator^, 16, ||, 

S A (q 2 ) = exp (-^) . (8) 



Inserting eqs.(||, |8|) into eq.((^), we obtain 



T\ ' T2 

V(r) = a w — - — [ai ■ <T2VA{r,m w ) + Su{r)tA(r,m n )] 



+ -k-C 2 V 2 -^IV#-#k(r) 
M I r dr r dr J AV ' 

+ Z 1 Z 2 av A (r,0) (9) 
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with S±2(r) = 2>a\ ■ f a 2 • r — a\ • a 2 , 



^i(r) = /^e^ r S A (q 2 ), (10) 



Mr,m) = f^e^4^, (11) 



S\(—m 2 ) J (2tt) 3 q 2 + m 2 

r did 
m 2 dr r dr 



r did 

t A (r,m) = — —VA(r,m) (12) 



and 



,2 ^2 



a^S A (-ml)^. (13) 

The factor 5^" 1 (— m 2 ) in the definition of v A (r, to) is introduced so that v A (r, m) and £a(^ m) 
become Yukawa functions in the large-r limit: 

e -mr 

lim v A (r,m) = , 

Ar>l r 

/ 3 3 \ e -mr 

lim t A (r, m) = 1 + + 7 ^ . (14) 

Ar>i \ mr \ynry ) r 

In eq.(||), we have absorbed the 5\(r) term appearing in the one-pion-exchange potential 
into the nucleon contact interactions by redefining Cq as 

C = C - 2 , jzn ■ t 2 01 ■ a 2 . (15) 

Hereafter we will omit the prime on Co- As for the one-pion-exchange potential, we set 
a n = 0.075 following 0, 



Let us write the potential for each channel of interest. 

• The np and pp 1 So channel: Here tensor interactions are ineffective and so the 
potential reduces, with eq.(|^), to 

Vir^So) = -a^ A (r,m.) + ^[c 1S0 -A 2 (AV-3)C 2 1S0 ]5|(r) 

+ ZiZ 2 aw A (r,0). (16) 

The coefficients Cq |° differ in the pp channel from those in the np channel due to the 

charge-independence breaking interactions. 

#5 Note that our EFT potential resembles the potential used in Nijmegen group's multi-energy analysis 
[ ]17| . They used - in fitting their potential - 1787 pp and 2514 np accurate scattering data below Ti a b = 350 
MeV published between 1955 and 1992. They separated their potential into the long-range potential Vl and 
the short-range potential Vs such that V(r) = Vl(t) for r > b, and V(r) = Vs{r) otherwise, with b — 1.4 
fm. Their long-range potential consists of electromagnetic interactions and one-pion-exchange potentials, as 
well as heavy-meson-exchange (such as uj, p and rj) potentials, but since Vl is applied only for r > b, the 
heavy-meson exchange is irrelevant except for achieving high precision in the large momentum region. The 
short-ranged part is parameterized by a square-well potential (of radius b) with an ^-dependent depth. The 
Independence is given by a finite series in increasing powers of E. Taking 1^3 terms for each channel 
(except for the pp 1 So channel for which four terms are retained), they achieve a very accurate fit with 
X 2 /data = 1.08. 
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The deuteron 3 Si — 3 Di coupled channel: The potential is of the form 



F(r; 3 Si- 3 £>i; 



V${r) + S 12 {f)V${r) 



(17) 



where, due to eq.(|l 

V$(r) 
VS(r) 



4tt 



-a^ v A (r, m n ) + — [Cf 1 - A^(AV - 3) Cf l \ <ft(r), 
-a n t A (r, m7T ) - ^^=- AV(5i(r). 



(18) 



With the regulator (||), it is a simple matter to work out the explicit forms of the 
quantities we need: 

d 1 d 



V^i(r) = A 2 (AV - 3) Sl(r), r---6l(r) = AV S 3 A (r) 



and 



<t(r) 
v A (r, m) 
v A (r,0) 



A 3 



exp 



(2tt)I 



1 

27 



-Ar + ^ 



(r 



where erf(rc) and erfc(x) = 1 — erf(x) are the error functions, 

erf(x) = —i= [ dte~ t . 
yvr Jo 



(19) 

(20) 

(21) 
(22) 

(23) 



Finally the radial functions, uo(r) for 1 Sq and u d (r) and Wd{r) for 3 S\ — 3 -Di, are 
obtained by solving the Schrodinger equations 



and 



4,'(r) + M - V(r; l S )} u (r) = 0, (24) 



v%(r)+M[E-Vg(rj\ud(r) = V8MV^(r)w d (r), 
w'^r)-^w d (r)+M\E-Vg(r) + 2V£(r)]w d (r) = V8MVS(r)u d (r), (25) 



with the boundary condition 



«o(0) = u d (0) = ta d (0) = 0. 



(26) 



One caveat is in order here. Solving the Schrodinger equation with the irreducible 
vertex calculated to NLO is consistent with the Weinberg scheme. As pointed out in 16, 0, 
however, it is not strictly consistent with the chiral counting order by order since the 
Gaussian cut-off function brings in all powers of p/A ~ Q/A and so does the iteration to 
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all orders of the irreducible vertex that is involved in the evaluation of reducible diagrams. 
Nevertheless, we believe that our procedure is safe from uncontrollable errors. First of all, 
in our treatment, the cutoff has a physical meaning and should not be sent to infinity. 



Therefore there is no problem of the sort described by Beane et a/.[16| and by Luke and 



Manohar [18], which one would encounter if one were to send the cutoff to infinity or if one 
were to use the naive dimensional regularization (DR) in the Weinberg scheme. Furthermore 
to the order we work in computing the irreducible diagrams, our method is consistent with 
the power counting. Possible errors reside only in higher-order terms that are being dropped. 
A similar conclusion was reached by Gegelia J2T 

3 Renormalization 

For a finite cutoff A, only finite quantities enter. However renormalization is required 
to render the theory consistent (see Q for details). In particular, certain cutoff-dependent 
terms have to be absorbed into the constants before determining them from empirical data. 
We proceed as follows. For a given cutoff and for each channel, we determine the constants 
C's and D's in (|5|) by relating them to the low-energy scattering parameters appearing in 
the effective-range formula. 

• The np 1 Sq channel: The effective-range formula for this channel is 

p Cot 5 = -— + lr p 2 + O(p 4 ) (27) 
«o 2 



where p = V ME is the relative momentum and the subscript '0' denotes the np 1 Sq channel, 
i.e., 5q = 1 5 np and similarly for the scattering length ao and the effective range r^. The 
low-energy scattering parameters, ao and ro, can be given by the radial wavefunction uo(r) 
at zero energy: 

ru' (r) — u$(r 



«o 



(28) 

r>R 



n' (r) 

r = 2 / dr <Pl(r)-u^r) , (29) 



oo 



where R denotes the range of the potential such that V(r)=0 for r > R, and (f>o(r) is the 
large-r asymptotic wavefunction of uo(r), 

T 

lim uo(r) = 4>o{r) = 1 . (30) 

r-*oo ao 
The experimental values of ao and tq which we shall use to fix Cq P and Crj 1 are: 

a^P = -23.749(8) fm, r[j xp = 2.81(5) fm. (31) 

• The pp 1 Sq channel: In this channel, the infinite-ranged Coulomb potential must 
be taken into account. There are some subtleties in doing this because there are many 



9 



different definitions of the low-energy scattering parameters depending on how the long- 
ranged 0(a 2 ) corrections are treated [19]. One reliable scheme was described in full detail 



in ]pT| 1, which we will follow here. We define the phase shift 6 - and consequently the 
low-energy scattering parameters a c and r - in the absence of the potentials of order 
0(a 2 ) and take into account these higher order EM effects by perturbation theory. The 
expressions obtained in |ll]] for the low-energy scattering parameters are 

a 1 "c+nW^W-uc+nW^W 



a 



(32) 

r>R 



r c = 2 / <Ir 

with 



m p a it' c+N (r)0i(r) -mc+nW^'iW 

oo r 

0cW-«c+nW (33) 







Mr)=Mr)-^^c, (34) 
m p a a° 

where m p ~ M is the proton mass, uc+i<i(r) is the radial wavefunction that approaches the 
function 4>c( r ) asymptotically, and <j>i(r) and 4 > 2{ r ) are the radial wavefunctions calculated 
with only the Coulomb interaction at zero energy, 

4>i (r) = 2^/m p arKi (2^Jm p ar) , 

</> 2 (r) = ^/m p arli(2^/m p ar). (35) 

Here -fta(-z) and ii(z) are the modified Bessel functions. In fixing Cq P and , we shall use 
the following experimental values (extracted from the Nijmegen multi-energy analysis Jl7[ ) 
of a c and r : 

a^ p = -7.8196(26) fm, rg p = 2.790(14) fm. (36) 

• The deuteron 3 Si — 3 D\ channel: In this channel, we have the choice of renor- 
malizing the theory either at zero energy or at the deuteron pole position E = —Bd, where 
Bd is the deuteron binding energy. Here we shall adopt the latter. The coefficient D% will 
be determined by the experimentally well-known deuteron D/S ratio rjd- Near the deuteron 
pole position, the effective-range formula has the form 

1 

2' 

where <5<j is the 5-wave eigenphase and 7 = yJMBd is the deuteron wave number. Thus the 
quantities Bd-, Pd and rjd are inputs to determine the three constants Cg , and D^. Since 
by construction our potential is energy-independent, the wavefunction normalization factor 
A s is not an independent quantity. It is related to the input parameters by 



P cot 5 d = -7 + ^p d (p 2 +l 2 ) + (p 2 + 7 T (37) 



„2 , _,2\2 



A '"(1 + ^)(1-7W)- (38) 



In our numerical calculation, we shall use the following experimental values for Bd, A s and 
Vd- 

B° xp = 2.224575(9) MeV, Af p = 0.8846(8) fnJ, rj e d xp = 0.0256(4). (39) 
From eq.pl) we get pd = 1.7635(46) fm. 
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4 Results 



As stated, we shall use the effective-range parameters to fix the coefficients of (|5|) 
for the np and pp 1 So channels and Bj, A s and rjd for the deuteron channel. Given solu- 
tions for the Schrodinger equation with the nucleon-nucleon potential thus obtained and 
the electroweak current derived to the leading order, we are in a position to carry out a 
parameter-free calculation of the Ml matrix element for the np capture process and the 
Gamow- Teller matrix element for the pp fusion process. Furthermore, the deuteron static 
properties such as the charge radius r^, the quadrupole moment Qd and the magnetic mo- 
ment /j,d can be predicted. The precise value of the cutoff is of course not determined by 
the theory but the consistency of the scheme requires that the results should be insensitive 
to its precise value provided that it is picked judiciously: e.g., around the pion mass if the 
pion is integrated out and around the two-pion mass if the pion is present. 

The relevant electroweak matrix elements we shall focus on are 



oo 



Mmi = / dru d (r)u (r), (40) 
Jo 

poo 

M G t = / dru d (r)u c +N(r). (41) 



o 



The 0(a 2 ) (VP and C2) corrections to the GT matrix element are calculated perturbatively |To|l 
using the potentials 



V§ 2 (r) + V^ 2 (r) = -^L(^ + V N (r)), (42) 



r 



where m e is the electron mass, and Vv( r ) denotes the potential ( |lq ) but without the 
Coulomb interaction. 

4.1 Without the pion 

We first discuss the case where the one-pion-exchange term in eq.(|5j) is absent. The 
results with the NLO terms are given in Table |l|. There the scattering lengths and effective 
ranges for the np and pp channels and Bd, A s and r\d for the deuteron channel are used as 
inputs to determine the constants of (||). We see that the results as a whole are stable against 
changes in the cutoff as long as the cutoff is taken above the pion mass. The theoretical 
values are compared with experimental data when available and with the Argonne v\% 
results^ otherwise. There is good agreement between the calculated and experimental 
values. 



* 6 We define these matrix elements with all other factors stripped off. In |H and juj, .Mmi an d A^gt were 
denoted by Mib and , respectively. 

* 7 As we discuss below, the results given by the Argonne Vis wavefunctions are fit to experiments and 
hence can be considered as "experiments" . 



11 



Table 1: The next-to- leading order (NLO) results without the pion field. The static properties of 
the deuteron, the Ml and GT amplitudes and the VP and C2 contributions to the GT amplitude 
are listed for various choices of the cutoff A. The low-energy input parameters are the scattering 
lengths and effective ranges for the np and pp So channels, and B^, A s and rjd for the deuteron 
channel. Their experimental values are given by eqs.(|3l], |3(], |39| ). 



A (MeV) 


140 150 175 200 225 250 


Exn. Uisll 1 


r d (fm) 
Qd (fm 2 ) 
Pd (%) 


1.973 1.972 1.974 1.978 1.983 1.987 
0.259 0.268 0.287 0.302 0.312 0.319 
2.32 2.83 4.34 6.14 8.09 9.90 
0.867 0.864 0.855 0.845 0.834 0.823 


1.966(7) 1.967 
0.286 0.270 
5.76 

0.8574 0.847 


Mmi (fm) 


3.995 3.989 3.973 3.955 3.936 3.918 


3.979 


Mqt (fm) 

*Xf (%) 
^ B 2:C (%) 
^B 2:N (%) 


4.887 4.881 4.864 4.846 4.827 4.810 
-0.45 -0.45 -0.45 -0.45 -0.45 -0.44 
0.03 0.03 0.03 0.03 0.03 0.03 
-0.19 -0.19 -0.18 -0.15 -0.12 -0.10 


4.859 
-0.45 
0.03 
-0.21 



These results incorporate up to NLO terms. Now, what happens when one limits the 
calculation to the leading order? The answer is that it does not work. It is not difficult to 
understand why this is so. Doing a leading order (LO) calculation corresponds to ignoring 
the q 2 term in (||) and, without the q 2 term - and without the pion, there will be no way 
to reproduce all the deuteron properties as well as the electroweak matrix elements. As 
shown in Fig. [I] for fixed CM momenta p ~ and p ~ m w , the low-energy np 1 Sq 

phase shifts calculated in LO (denoted by the dot-dot-dashed curve in Fig. |]) depend 
strongly on the cutoff. This strong cutoff dependence, however, is considerably reduced, 
when the NLO term is included (dashed curve). The improvement is more dramatic, the 
lower the momentum probed is. The introduction of the pion field makes the range of 
the relevant cutoff considerably greater, approaching the next momentum scale involved, 
namely, > 2m n . Thus the consistency condition for a viable effective theory with a given 
set of effective degrees of freedom requires going to sufficiently high order - NLO in the 
present case. To what order the calculation needs to be pushed and what additional degrees 
of freedom are required will of course depend on the problem one is addressing. 

4.2 With the pion 

We shall consider the cases with and without the 0(q 2 ) term in (0), the former 
corresponding to LO and the latter to NLO. 

• The LO calculation: Here the scattering lengths for the 1 5'o channel and the 
deuteron binding energy are inputs. The resulting predictions are given in Table [2|. 
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Figure 1: The np 1 Sq phase shift (degrees) vs. the cutoff A for fixed CM momenta p = 70 MeV 
(left) and p = 140 MeV (right). The next-to-leading order (NLO) results are given by the solid 
(with pions) and dashed (without pions) curves, and the leading-order (LO) results by the dot- 
dashed (with pion) and dot-dot-dashed (without pion) curves. The horizontal line represents the 
experimental data obtained from the Nijmcgcn multi-energy analysis. The NLO result without pions 
is drawn only up to A > 300 MeV, above which the theory becomes meaningless, that is, unnatural 
as well as inconsistent. 



With the pion, the deuteron does exist with its properties qualitatively reproduced. This 
is of course not surprising to nuclear physicists since the significant role of the pion in the 
deuteron structure has been known for a long time. What is surprising, however, is that 
even the electroweak matrix elements are not far from the "experimental values", though 
some 4 ~ 5 % higher. The effective ranges are not controlled without the NLO potential, 
so it is clear that the phase shifts can be reproduced only for very low energies. Actually 
we observe in Fig. |l| that while the LO result with the pion (denoted by dot-dashed) is 
better than the LO results without the pion (dot-dot-dashed curve), it is much worse than 
the NLO result without the pion (dashed curve). This feature, visible also in Fig. |3|, will 
be discussed in more detail later. An important point to note here is that when the pion 
field is present, the cutoff climbs up to roughly twice the pion mass. This is consistent with 
our expectation. 

• The NLO calculation: The pion contribution is fixed by chiral symmetry, so the 
parameters of the theory are determined in the same way as for the NLO calculation in the 
absence of the pion. The results are given in Table g 

The presence of the pion in the NLO calculation brings a noticeable improvement in 
the accuracy of the prediction making the dependence on the cutoff markedly weaker. This 
is true not only for the static properties (Table ^) but also for the phase shift for p ~ 
(see the right panel of Fig.[l|). The cutoff required is roughly twice that of the pionless 
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Table 2: The leading order (LO) results with pions. Note that only the scattering lengths for the 1 S'o 
channel and the deuteron binding energy are inputs. All the quantities except Co's are predicted; 
in particular, the wavefunction normalization A s and the D/S ratio rjd are not input but predicted 
here. 



A (MeV) 


300 400 500 600 700 


Exp. vis 11] 


C (ImJ 
C P P (fm) 
O (imj 


— U.ozo — U.ozl — U.oUl — U.4zz — U.oDD 
-0.810 -0.615 -0.500 -0.423 -0.367 
n Q7Q n n i /m n 70^ n 7sq 




A s (fin - 3) 

m 

r d (fm) 
Qd (fm 2 ) 
Pd (%) 


0.853 0.854 0.861 0.866 0.869 
0.0227 0.0250 0.0260 0.0263 0.0265 
1.897 1.902 1.918 1.930 1.936 
0.218 0.249 0.266 0.274 0.277 
2.65 4.59 5.97 6.73 7.08 
0.865 0.854 0.846 0.841 0.839 


0.8846(8) 0.885 
0.0256(4) 0.0250 
1.966(7) 1.967 
0.286 0.270 
5.76 

0.8574 0.847 


Mux (fm) 
Mgt (fm) 


4.176 4.186 4.164 4.140 4.122 
5.054 5.063 5.042 5.019 5.002 


3.979 
4.859 



case and this again is consistent with our expectation. Note that the charge radius, the 
quadrupole moment and the magnetic moment of the deuteron all come out very close to 
the experimental values, indicating that exchange-current corrections, if any, should be tiny. 
As remarked in the discussion section, this is expected on a general ground. The Ml matrix 
element and the GT matrix element receive a few per cents exchange-current corrections 
which are well controlled for the former and somewhat less well controlled for the latter. 
See later on this matter. 

To give a more quantitative idea as to how much the incorporation of the pionic 
degree of freedom reduces the cut-off dependence of our calculation, it is informative to 
consider £ (Ml) and £ (GT) defined by 



mil.^ ^b^ , £(GT)^ M ^ (43) 

JVI M1 JVI GT 



where -M]^ and denote, respectively, the Ml transition matrix element of our NLO 

calculation and that of the Argonne potential, and similarly for £(GT). Here we are taking 
the Argonne potential as "experiment" (see the next subsection for more explanation on 
this point). In Fig. || are plotted the deviations from the "experiment" as a function of 
the cutoff. We see that as in the case of the phase shift for p ~ m n , the inclusion of the 
pion degree of freedom markedly reduces the cutoff dependence, in conformity with the 
requirement for a viable effective theory that the more accurate the power expansion, the 
less cutoff-dependent physical observables should become. 
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Table 3: The NLO results with pion field. See the caption of Table ffl. 



A (MeV) 


200 250 300 350 400 500 


Exp. u l8 [ll 


fd (fm) 
Qd (fm 2 ) 
Pd (%) 


1.963 1.965 1.966 1.967 1.968 1.968 
0.261 0.268 0.272 0.273 0.274 0.274 
3.16 4.11 4.77 5.16 5.35 5.39 
0.862 0.856 0.853 0.850 0.849 0.849 


1.966(7) 1.967 
0.286 0.270 
5.76 

0.857 0.847 


A^mi (fm) 


3.987 3.976 3.968 3.963 3.958 3.952 


3.979 


Mgi (fm) 

^B P (%) 
^B 2:C (%) 
^B 2:N (%) 


4.884 4.874 4.867 4.862 4.859 4.854 
-0.46 -0.45 -0.45 -0.45 -0.45 -0.45 
0.04 0.03 0.03 0.03 0.03 0.03 
-0.24 -0.20 -0.15 -0.14 -0.21 -0.42 


4.859 
-0.45 
0.03 
-0.21 




-1%- 

1 ^ 

100 200 300 

A (MeV) 

Figure 2: £(M1) (upper) and £(GT) (lower) vs. the cutoff A 
results with pions and the dotted curves without pions. 



400 500 



. The solid curves represents the NLO 
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4.3 Comparison with a realistic phenomenological potential 

In our earlier calculations of the np capture process ||] and the proton fusion process 



[10], we resorted to the wavefunctions obtained from the Argonne vis potential - which 
is one of the most realistic phenomenological potentials fit to all available two-nucleon 
data up to TL a b=350 MeV - to calculate the matrix elements of the Ml and GT operators 
derived to the next-to-next-to leading order (NNLO) in chiral expansion. It might appear to 
some purists that this "hybrid" approach bringing in information from outside of the chiral 
perturbation theory framework is not consistent with the premise of an effective theory 
and hence throws doubt on the value and validity of such calculations. In this subsection, 
we argue that such a doubt can be dispelled because the hybrid model turns out to be in 
precise numerical agreement with the effective theory scheme. For our purpose, it suffices 
to limit our consideration to the leading-order matrix elements. 

In H], we showed that the leading-order Ml matrix element calculated in the hybrid 
method was is in good agreement with the result of the effective field theory approach in 
the absence of the pion. Here we confirm that the Ml matrix element as well as the GT 
matrix element (figuring in the pp fusion process) calculated both with and without pions 
quantitatively support the equivalence of the hybrid approach and the effective field theory 
approach with cutoff. We see from the NLO calculations given in Tables [l] and || that the 
electroweak matrix elements calculated in NLO in the effective field theory agree with those 
computed with the Argonne vis within 1 % error, with very little dependence on the cutoff. 
With the pion included the agreement is even more striking as one can see in Figure |2[ 

Comparison is made also for the 0(a 2 ) EM corrections to the pp fusion matrix ele- 
ment. These are obviously small corrections, smaller than higher-order power corrections 
ignored in the calculation but they cannot be ignored in the solar proton burning process as 
discussed in p0]| . As one can see in Tables [[] and |, the agreement is again almost perfect. 

5 Meson-Exchange Currents and Short-Range Correlation 

So far we have considered the matrix elements of only one-body electroweak cur- 
rents. In considering meson-exchange currents that appear as higher-order corrections, it 



is useful to recall the "chiral filter mechanism" first studied in [15, 21, 22]. The chiral 
filter mechanism is the statement that, whenever soft-pion exchange is not forbidden by 
symmetry or kinematically unsuppressed, the dominant term is given by chiral symmetry 
and the corrections to it should be suppressed. A corollary to this is that if soft-pion ex- 
change is forbidden by symmetry or kinematically suppressed, then higher-order terms enter 
non-negligibly so that the predictive power of the standard chiral perturbation strategy is 
appreciably weakened if not lost. 

According to this argument, we can make the following qualitative statements: 
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• The corrections to the single-particle Ml matrix element in the np capture process™ 
are given by one-soft-pion exchange (NLO) and loop corrections (NNLO). Since the 
LO one-body operator comes from the vertex with z/j = 1, these corrections are, re- 
spectively, 0(Q) and C(Q 3 ) relative to the LO operator. Working to 0(Q 3 ) relative 
to the (non-vanishing) leading order, we need to have wavefunctions (or potentials) 
accurate up to C(Q 3 ) for the one-body matrix element, to 0(Q 2 ) for the soft-pion 
exchange and to O(Q ) for one-loop correction. Since parity allows only even powers 
in the potential, it is sufficient therefore to have an 0(Q 2 ) accuracy. Our effective 
field theory to NLO meets this requirement. Let A^j^ 3,1 = A1mi(1 + Tty be the ma- 
trix element of the full Ml operator, with 1Z being the ratio of the matrix element 
of the exchange current operator to that of the leading Ml operator, AImi- With a 
"low-quality" phenomenological potential, we will very likely get a poor estimate of 
the leading one-body matrix element A^mi mostly because of its inability to repro- 
duce the low-energy scattering parameters. Even so, as explained above, the ratio 
1Z calculated with the same wavefunction is expected to be much less sensitive to 
the quality of the phenomenological potentials. That is, 1Z should be to a large ex- 
tent model- independent. Another way of arriving at this conclusion is as follows. As 
we have shown in this article, the one-body matrix element A^mi can be calculated 
with great accuracy with the NLO EFT potential as well as with a "high-quality" 
phenomenological potential such as the Argonne vi$, both of which reproduce all the 
relevant low-energy scattering observables with high precision. Since the total matrix 
element Af]yQ al should be the same regardless of whether we adopt the EFT frame- 
work developed here or the Argonne vis wavefunctions, we conclude that the ratio 1Z 
should not depend on which of the two methods is used. 

One non-trivial consequence of this argument is that the ratio 1Z should not depend 
upon how the current is regularized. Since the NLO correction term is made of tree 
diagrams only, the ratio should be reliably given by the Argonne v±s wavefunction. As 
for the NNLO term, it was found in || that all but one term (hereafter referred to as k 
term) is fixed by experiments. This n term which is left undetermined by experiment 
arises from a NNLO two-body counter term in the current. If one were to use the 
hybrid method of || , this k term - which is of zero range - could not contribute since 
it is suppressed by the short-range correlation in the wavefunction. Furthermore all 
zero-range operators that would appear as higher-order counter terms would also be 
suppressed. On the other hand, if one were to naively calculate the matrix element 
of the same zero-range operators in an EFT with the cut-off regularization used in 
this paper, the contribution would be non-zero unless the constants (e.g., k) happen 
to be all zero. So it would seem that the two approaches give different results for 

# 8 The same discussion applies to the time component of the axial current, that is, the axial-charge tran- 
sitions (first-forbidden /3-decays) in nuclei. For detailed discussions, see p3|. 
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non-zero re's at NNLO and higher orders. We conjecture that this difference can be at 
least partially resolved by appealing to the operator product expansion (OPE) in the 
wavefunction discussed in the context of EFT by Lepage 0. The OPE implements 
in the wavefunction the physics of perturbative QCD which is left out in the EFT 
framework discussed here. For long-wavelength operators (such as the single-particle 
Ml and GT operators), the effect of OPE would not be important. However, operators 
probing short-distance physics such as the zero-range k term should be affected by the 
OPE factor in the wavefunction. Our conjecture is that the short-range correlation in 
the standard nuclear physics approach and the OPE in the EFT approach a la Lepage 
capture qualitatively the same physics. At present we do not know how to make the 
connection between the two descriptions. We suggest that the lack of knowledge on 
the precise connection accounts for the uncertainty in the former approach manifested 
in different results for different forms of short-range correlation functions and for 
different short-range correlation sizes. This uncertainty - small in some cases like the 
np capture - is the price to pay in "killing" the higher-order counter terms, which we 
expect to be unimportant for very low-energy processes we have been considering but 
significant for higher-energy probes. 

For completeness we mention that an NNLO calculation for the exchange current 
contribution performed using the Argonne v\8 wavefunctions gave ~ 10% contribution 
to the np capture cross section leading to a precise agreement with the experiment 
The major contribution was given by the dominant soft-pion term, with only about 1 
% coming from NNLO loop corrections. The uncertainty in the calculation was about 
1 % in the cross section, representing the above-mentioned uncertainty in the short- 
range correlation that was used to "kill" the undetermined delta-function counter 
term. 

The power counting of the GT matrix element is slightly different from that of the 
Ml operator. Here the one-body GT operator is not suppressed by symmetry but 
one-soft-pion-exchange current is. As a consequence, the leading correction (one- 
pion-exchange) is of order of C(Q 3 ) while loops enter at 0(Q 4 ) or higher compared 
to the one-body operator. This means that the tree-order contributions involving one 
pion-exchange two-body operators and loop contributions and/or higher-dimensional 
counter terms that come at next order must be comparable in importance. So far 
a consistent calculation exists only in tree order, which predicts a ~ 4% correction 
(an enhancement) to the leading GT matrix element squared|l(|. The (corollary to 
the) chiral filter mechanism implies that this correction could receive a substantial 
contribution from higher-order or short-distance corrections. We should mention that 
the Landau- Migdal correction to g& in medium discovered a long time ago| 23||^ 
corresponds to this class of effects. 



A more modern interpretation in heavy-baryon chiral perturbation theory is given in 
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The exchange-current corrections to the deuteron properties, e.g., the charge radius, 
the quadrupole moment and the magnetic moment, are also suppressed by the chiral 



filter [15,^]. An additional suppression occurs due to isospin symmetry. This explains 
why in Tables |l| and ||| the theory without exchange current corrections agrees - almost 
precisely - with experiments without any further corrections. 



6 Breakdown of Effective Field Theory? 

Given a cutoff and a set of degrees of freedom chosen to work with, where is the given 
theory expected to break down and become inaccurate? By way of answering this question, 
let us look at low-energy np 1 Sq phase-shifts. In Fig. || are plotted the phase shifts vs. the 
CM momentum p for the lower and upper limits of the cutoff both in the absence and in 
the presence of the pion degree of freedom. The range A = (100 ~ 300) MeV represents the 
range of the cutoff in which the NLO theory without pions should be valid and the range 
A = (200 ~ 500) MeV is the corresponding range when the pion field is incorporated. 
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Figure 3: The np 1 Sq phase shift (degrees) vs. the center-of-mass (CM) momentum p. We show 
the predicted results both without (left) and with (right) the pion field. For each case, our NLO 
results are given for two extreme values of the cutoff: The solid curve represents the lower limit, 
A = 100 MeV (without the pion) and 200 MeV (with the pion), and the dotted curve the upper 
limit, A = 300 MeV (without the pion) and 500 MeV (with the pion). The "experimental points" 
obtained from the Nijmegen multi-energy analysis are given by the solid circles. Our theory with 
A = 150 MeV (without the pion) and A = 250 (with the pion) are shown by the dashed line to show 
that the theory is in almost perfect agreement with the data up to p < 200 MeV. 



There are several lessons one can learn from Fig. |||. By fiat, all the constants of 
the theory are fixed by the scattering length and the effective range, that is, the first two 
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terms in the effective-range expansion are fixed. However the Gaussian cutoff introduces 
terms of 0[(p/A) n ] for n > 3. The theory works fairly independently of the cutoff up to 
p ~ 100 MeV if the pion field is included, and up to p ~ 50 MeV if the pion is absent; 
beyond these limits the cutoff dependence becomes noticeable for both cases. Since the 
two-term effective-range formula fits data up to p ~ 140 MeV, higher-order terms must 
be responsible for the deviation from the experiment. Thus, in order to recover the fit, 
terms of order 0[(p/A) n ] with n > 3 must be included in the (irreducible) potential V. 
In fact the two-term effective-range formula works so well due to the fact that the third 



term of the effective-range formula, effective volume term, is very close to zero [25]. It then 
means that the two-term effective-range formula is essentially the same with the three-term 
effective-range formula with the effective volume term fine-tuned to zero. In order for this 
to happen in our theory, the contribution from the Gaussian cutoff and the corresponding 
higher-order counter term - which must be added - should cancel to give a vanishingly small 
volume term. Therefore it would be "unnatural" if our NLO calculation, without involving 
a "fine-tuning", were to work as well as the "three-term" effective-range formula. The left 
panel of Fig. 3 (NLO without the pion field) is thus "natural" . On the other hand, the right 
panel of Fig. 3 (NLO with the pion field) has the accuracy comparable to the "three-term" 
effective-range formula, which should be understood as the role of the pion field. 

Given that we are truncating the series, thereby making the theory manifestly ap- 
proximate, it seems reasonable to exploit the insensitivity to the precise value of the cutoff 
and do a "best fit." For this, we shall pick - within the range allowed - a specific value 
for cutoff in fitting experimental data. We should stress that this has nothing to do with a 
"fine-tuning" usually associated with an unnatural theory. From Table || we can see that, 
when the pion field is present, A ~ 250 MeV gives an overall fit of the deuteron properties 
as well as electroweak matrix elements. In Fig. |3| is plotted the phase shift prediction for 
this cutoff. The agreement with experiments is excellent up to p < 200 MeV- p^ Thus for 
this value of the cutoff, all low-energy two-nucleon properties are accurately described by 
the theory. 

We can do a similar "best-fitting" in the absence of the pion field and get a result 
which, though not as good as in the presence of the pion field, is still quite satisfactory as 
one can see in the left panel of Fig. ||| for A ~ 150 MeV. 



* 10 This accuracy of the effective-range formula to p ~ 140 MeV is somewhat unnatural from the point of 
view of EFT. It is perhaps fortuitous. Intuitively one would expect that the effective-range formula should 
work only up to p ~ m,/2, since the 5*-matrix element contains the left-hand cut starting at p 2 = — (m^/2) 2 . 

#11 In JtJ, an even more impressive agreement for the phase shift was obtained with the pion field included. 
However it should be noted that to achieve a better global fitting, the authors used values of the coefficients 
of the potential that are different from the ones determined by their renormalization procedure. As a 
consequence, the description for the threshold region became less reliable. 
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7 Discussion 



We have shown in this paper that the idea of effective field theory in its simplest 
form works remarkably well in nuclear physics of two-nucleon systems at low energy. We 
have used the cut-off regularization which renders the notion of effective theory physically 
transparent. It has been seen that the cut-off regularization also renders moot [Q] the issue 
of whether there is a breakdown in Weinberg's counting rule, which consists in applying 
the chiral counting only to irreducible diagrams, and subsequently summing up an infinite 
series of reducible graphs containing irreducible vertices of a specified order. 

We have further shown that the "hybrid approach" used in || to calculate the n+p — > 
d + 7 cross section in the framework of chiral perturbation theory (in the Weinberg version) 
is fully justified by a systematic effective field theory technique with a finite cutoff. 

One of the important by-products of the effective field theory for nuclei is precise 
calculations of nuclear weak processes that figure importantly in astrophysics. In [jnj , we 
computed in the "hybrid approach" the cross sections for the solar proton burning process 
p + p —> d + e + + v e to an accuracy of 1 %,^ T ^ and this has provided strong support for the 
estimate used by Bahcall and collaborators for the solar neutrino problem[2C, 26, p7|] . In 
particular it has resolved the possible controversy as to whether or not a field theoretic ap- 
proach would give a result basically different from the classic nuclear physics approach which 
is based on Schrodinger equation with realistic phenomenological N-N interactions f20(| . An- 
other potentially important by-product is that with the pion field incorporated, we should 
be able to treat meson-exchange currents to one-loop order (that is, 0(Q 4 ) relative to 
the leading GT operator). This would allow a model- independent calculation of processes 

likepl mi 

3 He + p -> 4 He + e + + v e (44) 

which may be measured at the SuperKamiokande and SNO. 

As often heralded by the practitioners of the dimensional regularization (DR), one of 
the advantages of DR over other regularizations is that the DR scheme preserves chiral in- 
variance and Lorentz invariance. Indeed, while giving an excellent result in a non-relativistic 
regime, the cut-off regularization as used here loses manifest chiral invariance and Lorentz 
invariance. This might make higher-order loop calculations difficult, if not impossible. 
Fortunately, in the tree-order treatment of meson exchange currents for the electroweak 
processes as well as the other two-nucleon processes considered in this paper, we did not 
have to deal with chirally non-invariant counter terms. However, this problem does show 
up at the level of loop calculations. We have not investigated this problem in detail but 



* 12 The main caveat in is in the meson-exchange contribution which comes out to be about 4 % when 
calculated to the chiral order C(Q 3 ). The higher order contribution starts from 0(Q 4 ), and may not be 
necessarily negligible compared with the 0(Q 3 ) contribution. So a more "conservative" error estimate may 
be 4 %, instead of 1 %. 

* 13 We would like to thank John Bahcall for calling our attention to this problem. 
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in calculating the process in (|44|), where meson-exchange current corrections are expected 
to be more important than in the cases we have considered so far, this issue will have to 
be resolved. A related issue was investigated by the authors of ||], who suggest that the 
difficulty is not insurmountable. Their work indicates a promising way of dealing with loop 



corrections in a finite cut-off field theory for the process (44) as well as for the proton fusion 
process. 
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Appendix 



We briefly comment on a recent discussion |7j on Weinberg's counting rule. For 



related discussions on this issue, see ||, |29|, [3C], [31, [3^, |33| . More recently, Birse, McGovern 
and Richardson |?4| pointed out that for two-nucleon scattering with a quasi-bound state 
nearby as involved in these discussions, fluctuations should be made around a non-trivial 
fixed point. We shall not address this issue here. 

The authors of Q have argued that the counting rule eq. (||) breaks down when applied 
to reducible graphs, which need to be summed to all orders to describe bound or quasi-bound 
states. If one were to use the usual dimensional regularization, this infinite series with a 
finite set of irreducible graphs for a potential would require an infinite number of counter 
terms to regularize the series 

In constructing counting rules of low-energy effective theories that are based on deriva- 
tive expansion, subtlety (if any) resides in loop integrals. To be definite, we consider an 
n-dimensional integral of a function /(p, q), where p stands for an external momentum 
and q a loop-momentum, 

/ A (p) = /Vq/(p,q), (45) 

where the superscript A indicates that the quantity is regulated with a cutoff A. Let /(p, q) 
be of order of Q u f when both p and q are of order of Q. Weinberg assumed that the integral 
should be counted as 

J A (p) = Q{Q n+v f). (46) 

Note that we do not write down here explicitly the power of the "large scale" A. If one 
naively assumes that I A (p) ~ A n Q u f ~ Q v f, where ~ stands for "the same order" , then one 
would find that Feynman graphs with arbitrary numbers of loops are not suppressed, which 
would make the effective theory completely useless. To elucidate Weinberg's assumption 
(|46|), let us expand I A (p) in decreasing powers of A, 

oo 

J A ( P ) = ^A«- fc / A (p) (47) 

k=0 



with 



4 a (p) ~ Q vs+k - (48) 



The first n-terms (0 < k < n) — which diverges when the cutoff A goes to infinity - 
appear as (finite) polynomial in p. If the loop integral is for an irreducible loop, then the 
polynomials are to be absorbed into the counter terms of the corresponding vertex functions 
in the renormalization procedure. Then these terms become "invisible" and hence do not 
affect any physical observables. This implies that, in constructing counting rules, one can 
effectively count the cutoff of the loop as order of Q, even though the "actual cutoff" A is 
much bigger than Q and is not directly related to Q. In our understanding, this feature 
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leads to Weinberg's prescription: Assign the "effective cutoff', A*, 

A* ~ Q < A, (49) 

to any loop integrals in irreducible diagrams. It is then to be expected that Weinberg's 
assumption ( filf ) cannot be simply extended to reducible loops. In fact, the assumption 
should not work for reducible loops if a system is to have a (quasi) bound state (see below) . 
Another problem is that, for reducible loops, it is not obvious how to absorb the polynomials, 
Ik(p) with k < n, into counter terms. The standard quantum field theory prescription tells 
us only how to renormalize vertex functions which receive contributions only from irreducible 
loops and counter terms. 

To further expound these difficulties, we examine the counting rules for reducible 
graphs. The Lippman-Schwinger equation is schematically of the form, T = V + VG°T, 
where V is a regulated potential, G° the free two-nucleon propagator and T the T-matrix. 
For a bound state, there is a pole in the T matrix satisfying the homogeneous equation 
T = VG°T, which implies VG° = 1. For quasi-bound states, we have VG° = 0(1). For 
simplicity, let us retain only the Co (leading-order contact) term and assume the separability 
of the potential. Then the required condition for a reducible loop graph is 

ci / ,f ' i4g '- a /*- 1 ' <5o) 

In finite cut-off effective theory, as noted in |J, Cq scales as A -1 and the integral / dq ~ A, 
so the condition ( |50| ) is trivially satisfied. In the Weinberg counting rule (|2|), as we have 
surveyed, every loop in n-dimensional space is counted as 0(Q n ). When reducible graphs 
are present, we need to introduce a modification to this rule. We note that the following 
rule is consistent with eq.(|5(]): Any reducible loop in n-dimensional space should be counted 
as 0(Q n ~ l ). It is also interesting to find that eq.([50|) implies that the very natural scale 
of | Co | ~ A" 1 forces any low-energy two-nucleon interaction to become non-perturbative 
provided that Co < 0, which is consistent with the shallow-binding of the deuteron and the 
almost-binding of np 1 Sq channel in nature. 

In the dimensional regularization (DR) scheme, however, realizing the condition ( |50"| ) 
is non-trivial due to the lemma J (f n q = 0. This could be a reason why EFTs based on the 
naive DR fails in describing low-energy ^o scattering @, H- In @, a new counting rule 



was suggested involving subtractions not only at D = 4 but for all D < 4, where D is the 
space-time dimension. In this scheme, the condition corresponding to (|50|) takes the form 



C /u ~ 1 (51) 

with 

Co = -rj~ ' ( 52 ) 

1/a — u 



where [i is the renormalization scale. Eq.(p2[) satisfies the corresponding consistency condi- 
tion eq.(|5l|), which is responsible for the success of the PDS scheme in contrast to the failure 
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of the naive dimmensional regularization.^ 15 However, if we were to count the scale [i as of 
order of Q as was done in jj]] , there would appear an inconsistency with the Georgi-Manohar 
counting rule 35] that ^jCo ~ -p ~ jp. On the other hand, if one were to adopt the PDS 
scheme but count the renormalization scale as [i ~ A appropriate for local operators, then 
one would be led to Co ~ ^ and hence to Weinberg's counting rule. 

* 14 One of the powers of the PDS scheme is the ability to write an analytic expression for Co in terms of 
the unnaturally large scattering length a, eq. (pS^) , reflecting a large anomalous dimensiion of the four- Fermi 
(baryon) interaction. This feature is somewhat buried in the cut-off method although it is implicit in it. 
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